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Abstract: We reconsider and update the QCD light-cone sum rules for B — > it form 
factors. The gluon radiative corrections to the twist- 2 and twist-3 terms in the correlation 
functions are calculated. The MS 6-quark mass is employed, instead of the one-loop 
pole mass used in the previous analyses. The light-cone sum rule for /^(g 2 ) is fitted to 
the measured (^-distribution in B — > nlui, fixing the input parameters with the largest 
uncertainty: the Gegenbauer moments of the pion distribution amplitude. For the B — > ir 
vector form factor at zero momentum transfer we predict /r_(0) = 0.26^g'^|. Combining 
it with the value of the product |V^b/^ (0)| extracted from experiment, we obtain \V u b\ = 
(3.5 ± 0.4 ± 0.2 ± 0.1) x 10~ 3 . In addition, the scalar and penguin B — > ir form factors 
IbM 2 ) and IbM 2 ) are calculated. 
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1. Introduction 

The form factors of heavy-to-light transitions at large energies of the final hadrons are 
among the most important applications of QCD light-cone sum rules (LCSR) JTJ. In this 
paper we concentrate on the B — > it transition form factors fg n , fg n and of the 
electroweak vector b — > u and penguin b — > (i currents, respectively. Previously, these form 
factors have been calculated from LCSR in (2|, ||, [I[ || ^, ^, |8|, ^, 10 1 , gradually improving 
the accuracy. 

The main advantage of LCSR is the possibility to perform calculations in full QCD, 
with a finite 6-quark mass. In the sum rule approach, the B — > tt matrix element is 
obtained from the correlation function of quark currents, rather than estimated directly 
from a certain factorization ansatz. This correlation function is conveniently "designed", 
so that, at large spacelike external momenta, the operator-product expansion (OPE) near 
the light-cone is applicable. Within OPE, the correlation function is factorized in a series 
of hard-scattering amplitudes convoluted with the pion light-cone distribution amplitudes 
(DA's) of growing twist. To obtain the B — > tt form factors from the correlation function, 
one makes use of the hadronic dispersion relation and quark-hadron duality in the S-meson 
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channel, following the general strategy of QCD sum rules []ll]] . More details can be found 
in the reviews on LCSR, e.g., in [12, 13, 14|. A modification of the method, involving 
.B-meson distribution amplitudes and dispersion relation in the pion channel was recently 
suggested in |ll| ; the analogous sum rules for B — > ir form factors in soft-collinear effective 
theory (SCET) were derived in pq[ . 

LCSR provide analytic expressions for the form factors, including both hard-scattering 
and soft (end-point) contributions. Because the method is based on a calculation in full 
QCD, combined with a rigorous hadronic dispersion relation, the uncertainties in the result- 
ing LCSR are identifiable and assessable. These uncertainties are caused by the truncation 
of the light-cone OPE, and by the limited accuracy of the universal input, such as the 
quark masses and parameters of the pion DA's. In addition, a sort of systematic uncer- 
tainty is brought by the quark-hadron duality approximation adopted for the contribution 
of excited hadronic states in the dispersion relation. Importantly, B — > tt form factors are 
calculable from LCSR in the region of small momentum transfer q 2 (large energy of the 
pion), not yet directly accessible to lattice QCD. 

The B — > irlvi decays, with continuously improving experimental data, provide nowa- 
days the most reliable exclusive V u \, determination. Along with the lattice QCD results, the 
form factor fg n {q 2 ) obtained [|io| from LCSR is used for the \ V u b\ extraction. Furthermore, 
the LCSR form factors f^(q 2 ) 

can provide inputs for various factorization approaches to 
exclusive B decays, such as QCD factorization (|T?j ], whereas the penguin form factor 
is necessary for the analysis of the rare B — > irl + l~ decay. Having in mind the importance 
of B — > tt form factors for the V u b determination and for the phenomenological analysis of 
various exclusive B decays, we decided to reanalyze and update the LCSR for these form 
factors. One of our motivations was to recalculate the 0(a s ) gluon radiative correction to 
the twist-3 part of the correlation function. Only a single calculation of this term exists 
H |i"0| , whereas the 0(a s ) corrections to the twist- 2 part have been independently obtained 
in Q and H. In what follows, we derive and present the explicit expressions for all 0(a s ) 
hard-scattering amplitudes and their imaginary parts for the twist-2 and twist-3 parts of 
the correlation function and some of these expressions are new. 

In the OPE of the correlation function the MS mass rn^jj) is used, a natural choice 
for a virtual 6-quark propagating in the hard-scattering amplitudes, calculated at large 
spacelike momentum scales ~ m^. Importantly, in the resulting sum rules we keep using 
the MS mass. Note that the value of m&(m{,) is rather accurately determined from the 
bottomonium sum rules. In previous analyses, the one-loop pole mass of the 6-quark was 
employed in LCSR. The main motivation was that the pole mass was used also in the 
two-point sum rule for the i?-meson decay constant /b, needed to extract the form factor 
from LCSR. In the meantime, the fs sum rule is available also in MS-scheme [18], and we 
apply this new version here. 

Furthermore, we fix the most uncertain input parameters, the effective threshold and 
simultaneously, the Gegenbauer moments of the pion twist-2 DA, by calculating the B- 
meson mass and the shape of f^ n {q 2 ) from LCSR and fitting these quantities to their 
measured values. In addition, the nonperturbative parameters of the twist-3, 4 pion DA's 
entering LCSR are updated, using the results of the recent analysis |ts[[. 
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The paper is organized as follows. In sect. 2 the correlation function is introduced and 
the leading-order (LO) terms of OPE are presented, including the contributions of the pion 
twist-2,3,4 two-particle DA's and twist-3,4 three-particle DA's. In sect. 3 the calculation 
of the 0(a s ) twist-2 and twist-3 parts of the correlation function is discussed. In sect. 4 
we present LCSR for all three B — > ir form factors. Sect. 5 contains the discussion of 
the numerical input and results, as well as the estimation of theoretical uncertainties, and 
finally, the determination of \V u b\. Sect. 6 is devoted to the concluding discussion. App. A 
contains the necessary formulae and input for the pion DA's. The bulky expressions for 
the 0(a s ) hard-scattering amplitudes and their imaginary parts are collected in App. B, 
and the sum rule for fs is given in App. C. 

2. Correlation function 

The vacuum-to-pion correlation function used to obtain the LCSR for the form factors of 
B — > ir transitions is defined as: 



F„ip,q) = ij d 4 x e l <i x {ir + ( P )\T{u(x)T^b(x),m b b(0)i l5 d(0)} |0) 
( F(q 2 , (p + q) 2 ) Pll + F(q 2 , (p + q) 2 )q^ , = 7^ 

I F T (q 2 , (p + q) 2 ) [p^q 2 - q^{qp)] , T„ = -ia^q 1 



(2.1) 



for the two different b — > u transition currents, For definiteness, we consider the Bd — > vr + 
flavour configuration and, for simplicity we use u instead of d in the penguin current, 
which does not make difference in the adopted isospin symmetry limit. Working in the 
chiral limit, we neglect the pion mass (p 2 = m 2 = 0) and the u-, (f-quark masses, whereas 
the ratio fjL n = m 2 /(m u + m^) remains finite. 

At q 2 <C m 2 and (p + q) 2 <C m 2 , that is, far from the 6-flavour thresholds, the b quark 
propagating in the correlation function is highly virtual and the distances near the light- 
cone x 2 = dominate. It is possible to prove the light-cone dominance, following the same 
line of arguments as in jl5|l . Contracting the 6-quark fields, one expands the vacuum-to 
pion matrix element in terms of the pion light-cone DA's of growing twist. The light-cone 
expansion [20] of the 6-quark propagator is used (see also ||): 



(01^(^(0)10) = -i 



d 4 k _ e _jfc. x . 



2 — k 2 



m 



1 

+g s J dvG^ a (vx) 



X a \ ij I jt + m 1 

I 2(m 2 - k 2 ) 2 a » v + m 2 -k 2 V%illv 



, (2.2) 

a/3 



where only the free propagator and the one-gluon term are retained. The latter term gives 
rise to the three-particle DA's in the OPE. Diagrammatically, the contributions of two- 
and three-particle DA's to the correlation function are depicted in Fig. [l|. In terms of 
perturbative QCD, these are LO (zeroth order in a s ) contributions. The Fock components 
of the pion with multiplicities larger than three, are neglected, as well as the twists higher 
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Figure 1: Diagrams representing the leading- order terms in the correlation function involving the 
two-particle (left) and three-particle (right) pion DA's shown by ovals. Solid, curly and wave lines 
represent quarks, gluons, and external currents, respectively. 

than 4. This truncation is justified by the fact that the twist-4 and three-particle corrections 
to LCSR obtained below turn out to be very small. 

In addition we include the 0(a s ) gluon radiative corrections to the dominant twist-2 
and twist-3 parts of the correlation function. The OPE result for the invariant amplitude 
F is then represented as a sum of LO and NLO parts: 



F(q\(p + q) 2 )=F (q 2 ,(p + q) 2 ) + 



a s C F 
4tt 



Fi(q 2 ,(p + q) 2 ), 



(2.3) 



and the same for F and F T . The leading-order (LO) invariant amplitudes Fq, Fq, and Fq 
including twist 2,3,4 contributions have been obtained earlier in |3|, ||, pi] ]. We present 
them here switching to the new notations of the twist-3, 4 DA's: 




, rnl + q 2 , ff , . m 2 ^(u) 



6m b [ ml-{ q + U pf\ Y ^ y 2 (m 2 b - (q + up) 2 ) 2 





(2.4) 
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Fo(q 2 ,(p + q) 2 ) = mbU 



du 



m i ~ (? + u vY 



' 6 



2 2 
m£ — q 

m? — (q + up) 2 



m b ~ (l + '"P) 2 



(2.5) 



W,(p + ,) 2 ) = ^ / ro g- ( t«p)^ (u) + 3K-T^) ^ (M) 







^ + ^2 7~ vo W( U ) 



2(m 2 — (g + up) 2 ) I 2 m 2 — (g + up) 2 / 

+m h U I dv f — — -2-2 < 2¥ 49r (o! i ) - (1 - 2t;)$ 47r (a 

J J \mt — (a + (at + aau 







[m 2 , - (g+ («i + a 3 v)p)' 

+2(l-2u)*4,r(ai)-$4*(c*)J, (2.6) 

where T>a = daidoL2da-z5{\ — a± — a<i — 03), and the definitions of the twist-2 ((p^), twist-3 
((frzir, 03j,-, $3tt) and twist-4 (^4^, tp4 n , $4^, ^4^, x J , 4 7r ) pion DA's and their parameters 
are presented in App. A. Note that all twist-4 terms are suppressed with respect to leading 
twist-2 terms, with an additional power of the denominator l/(m^ — {q + up) 2 ) compensated 
by the normalization parameter S 2 ~ A-qcd °^ the DA's. 

The calculation of the NLO amplitudes Fi,F\,Ff will be discussed in the next section. 

3. Gluon radiative corrections 



In the light-cone OPE of the correlation function (2.1) each twist component receives gluon 
radiative corrections. To obtain the desired NLO terms, one has to calculate the 0(a s ) 
one-loop diagrams shown in Fig. convoluting them with the twist-2 and two-particle 
twist-3 DA's, respectively. The diagrams are computed using the standard dimensional 
regularization and MS scheme. In addition, in our calculation the reduction method from 



1 22] is employed. 

The invariant amplitude -Fi in (|2.3j) is obtained in a factorized form of the convolutions: 



Fl{q 2 ,{p + q?) = f«J Q du^T^^p + q) 2 ^)^) 



T?(q 2 , (p + q) 2 , + T[(q 2 , (p + q) 2 ,u)^{u) , (3-1 
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Figure 2: Diagrams corresponding to the 0(a s ) gluon radiative corrections to the correlation 
function. 



where the hard-scattering amplitudes Ti, Tf' a result from the calculation of the diagrams 
in Fig|| The two other NLO amplitudes F\ and m^Ff have the same expressions with 

, respectively. The resulting expressions 
for all hard-scattering amplitudes are presented in App. B. Note that the LO expressions for 
the correlation functions in (2.4)-(p]q) also have a factorized, albeit a much simpler form, 
with the zeroth-order in a s hard-scattering amplitudes stemming from the free propagator 
of the virtual 6-quark. In particular, the twist-2 component in Fq is a convolution of 
T o = ml /[ml - (q + up) 2 ] with <p n (u). 

Let us mention some important features of the 0(a s ) terms of OPE. The currents u'y^b 
and m^bi^^d in the correlation function are physical and not renormalizable. Hence, the 
ultraviolet singularities appearing in T\ and T\ are canceled by the renormalization of the 
heavy quark mass. For an additional renormalization of the composite qc^b operator 
has to be taken into account. Furthermore, in the twist-2 term in (3.1) the convolution 
integral is convergent due to collinear factorization. As explicitly shown in [||, ||, the 
infrared-collinear divergences of the 0(a s ) diagrams are absorbed by the well known one- 
loop evolution [23] of the twist-2 pion DA. As a result of factorization, a residual dependence 
on the factorization scale fj,f enters the amplitude T\ and the twist-2 DA (p n . This scale 
effectively separates the long- and short (near the light-cone) distances in the correlation 
function. In the twist-3 part of Fi, the complete evolution kernel has to include the mixing 
of two- and three-particle DA's. To avoid these complications, and following ||, the twist-3 
pion DA's in (3.1) are taken in their asymptotic form: 4> p (u) = 1 and 4> a (u) = 6u(l — u), 
whereas the nonasymptotic effects in these DA's are only included in the LO part Fq. We 
checked that the infrared divergences appearing in the amplitudes Tf and T[ cancel in the 
sum of the (fi p and <j) a contributions with the one-loop renormalization of the parameter fj, n 
(i.e., of the quark condensate density). Finally, in accordance with j|, |l(]], all renormalized 
hard-scattering amplitudes are well behaved at the end-points u = 0, 1, regardless of the 
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form of the DA's. 

After completing the calculation of OPE terms with the LO (NLO) accuracy up to 
twist-4 (twist-3), we turn now to the derivation of the sum rules. 

4. LCSR for B — > ir form factors 



In the LCSR approach the B — > tt matrix elements are related to the correlation function 



(2.1) via hadronic dispersion relation in the channel of the bj^d current with the four- 
momentum squared (p + q) 2 . Inserting hadronic states between the currents in (|2.1| ) one 
isolates the ground-state l?-meson contributions in the dispersion relations for all three 
invariant amplitudes: 



F{q\(p + q?) 
F(q 2 ,(p + q) 2 ) 
F T (q 2 ,(p + q) 2 ) 



m 



B 



(P + 



+ ... 



m 2 B fB[fjM 2 ) + fBM 2 )\ 
m%-(p + q) 2 

2m 2 B f B flM 2 ) 
{m B + m 7r )(m^ 



+ ... 



+ • 



(4.1) 



where the ellipses indicate the contributions of heavier states (starting from B*ir). The 
three B — > it form factors entering the residues of the B pole in (|4.1[) are defined as: 



(7T + (p)\u %l b\B d (P + Q)) = VL(q 2 )P, + UbM 2 ) + fsM 2 )) ^ > 

iflM 2 ) 



{ir + (p)\ua tlu q u b\B d (p + q)) = q (2p^ + g M ) - (m B - roj^ 



mB + run 



(4.2) 
(4.3) 



and fs = (B(i\mbbi^d\0) /m 2 B is the S-meson decay constant. 

Substituting the OPE results for F, F and F T in l.h.s. of (4.1), one approximates 
the contributions of the heavier states in r.h.s. with the help of quark-hadron duality, 
introducing the effective threshold parameter s B . After the Borel transformation in the 
variable (p + q) 2 — » M 2 , the sum rules for all three B — > ir form factors are obtained. The 
LCSR for the vector form factor reads: 



f + B M 2 



e m%/M* 

2m|/B 



F (q*,M 2 ,s* 



+ ^L Fl {q 2 ,M 2 ,s B ) 



(4.4) 



where F ^(q , M , Sq) originates from the OPE result for the LO (NLO) invariant ampli- 
tude F 0(1) (q 2 ,(p + q) 2 ). 

The LO part of the LCSR has the following expression: 

.2 ,,2 „BX__2, / , J ...-^£/v w («) 



F (q 2 ,M 2 ,s%)=m£f n 



du e 



mb\ 6 L u 



m\ + q 2 



mt 



du 



+- 



m 2 u cP<fr&ir(u) 



mi 



4(m 



») du 2 



2 ( 

\m b f n 



u 



utp AlT (u) + J dv7p i7T (v) - h^u) \ >, (4.5) 

' ' 



-7- 



where u = 1 — u, uq = (m 2 — q 2 )/{sQ — q 2 ) and the short-hand notations introduced for 
the integrals over three-particle DA's are: 



(u 1 
I dal I 



dv ^3vr(aj) 
1 ' i'(-'H)/(l-ai) 

1 

2^ 4 7r(ai) - $47r(a;) 



«2 = 1 — Oil — 03, 

03 = (u — ai)/v 



Ii7T ^ = ~<h l ( j dai J V 



(u- ai )/(l- ai ) 

+2* 47 r(ai) - <*>47r(ai 



a-i = 1 — ol\ — 03, 
«3 = (u — ai)/-y 



The NLO term in ( |4.4[ ) is cast in the form of the dispersion relation: 
F l {q 2 ,M 2 ,s*) = ±- | dse-^Im^iq 2 ^) 

S B 

1 ( 

= — dse~ s / M / dul Im s Ti(g 2 , s, u) ^^{u) 
k J Jo I 



(4.6) 



Hit 



lm. s T'f(q ,s,u)4% w (u) + lm s T[(q ,s,«)^ 7r 



(4.7) 



where the bulky expressions for the imaginary parts of the amplitudes Ti,Tf,Tf are pre- 
sented in App. B. 

The LCSR following from the dispersion relation for the invariant amplitude F in ((OJ) 
reads: 



m%f B 



where 



F (q 2 ,M 2 ,s 3 )=mif 7T I due 



h{q\M 2 , S Z) + ^F l {q 2 ,M 2 ,s*) 



(4.8) 



mj, \ u 6u du 



'to 



+- 



2 _ „2 



^4ir(«) 



(4.9) 



Here the contributions of twist- 2 and of three-particle DA's vanish altogether. Combining 
(|4.4j) and (4.8) one is able to calculate the scalar B — > 7r form factor: 



/L(? 2 ) = /^(<z 2 ) + 



2 9 

— m~ 



(4.10) 
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Finally, the LCSR for the penguin form factor obtained from the third dispersion 
relation in ( |4.1| ) has the following expression: 



JbAQ ) — — 



where 



and 



2m%f B 



F^ q \M\ S Z) + ^Fnci\M\sE) 



1 

F^(q 2 ,M 2 ,s^) = m b f 7T J due 

f o 



1 ( 1 dfa n (u) _ m\u d 2 (j) i7T (u) _ t t ( \ \ \ 
m 2 b -q 2 \4 du 2(m 2 -q 2 ) du 2 47rW j J 



uM z < — 

u 3(m^ — q 2 ) du 



(4.11) 



(4.12) 



u 



du 



da 



dv 
v 



(«-ai)/(l-ai) 



2tf 4 ,r(ai) -(l-2u)$4* (a*) 



+2(l-2u)^(ai)-$4^(ai) 



«2 

a 3 



1 - ai - a 3 , 
(u — ai)/v 



(4.13) 



The NLO parts Fi and m^F^ in LCSR (4.8) and ( |4.11| ), respectively, are represented in 
the form similar to (POT), and the corresponding imaginary parts are collected in App. B. 



For fs entering LCSR we use the well known two-point sum rule 1 24] obtained from 
the correlator of two mbqi^b currents. The latest analyses of this sum rule can be found 
in |18|, 25]; here we employ the MS version [18]. For consistency with LCSR, the sum rule 
for Jb is taken with 0(a s ) accuracy. For convenience, this expression is written down in 
App. C. 

Note that the expressions for LCSR in LO are slightly modified as compared to the ones 
presented in the previous papers. We prefer not to use the so-called "surface terms" , which 
originate from the powers of l/(m 2 — (q + up)) n with n > 1 in the correlation functions. 
Instead, we use a completely equivalent but more compact form, with derivatives of DA's. 

The twist-2 NLO part of LCSR for , hence, the expressions for T\ and ImTi in 
App. B, after transition to the pole scheme (the additional expressions necessary for this 
transition are also presented in App. B) coincide with the ones obtained in pj. We have 
also checked an exact numerical coincidence with the twist-2 NLO part of the sum rule 
in ||, written in a different analytical form. The explicit expressions for the amplitudes 
Tf ,a , T\, Tf ,a , and T± \T^ p,a and their imaginary parts presented in App. B are new. The 
0(a s ) spectral density entering the LCSR for ft is given in 0] in a different form, that 
is, with the tt-integration performed, making an analytical comparison of our result with 
this expression very complicated. The numerical comparison is discussed below, in sect. 6. 
Furthermore, in [^] the LCSR for the form factor was obtained, and the imaginary 
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Figure 3: Replacing the integration intervals by the contours in the complex planes of u and s 
variables in the alternative procedure of the numerical integration of NLO amplitudes. 

part of T\ was presented. A comparison with our expression for ImTi reveals, however, 
some differences. 

Since the imaginary parts of the hard-scattering amplitudes have a very cumbersome 
analytical structure, we carried out a special check of these expressions. Each hard- 
scattering amplitude Xi, ... taken as a function of u, q 2 ,{p + q) 2 was numerically compared 
with its dispersion relation in the variable (p + q) 2 = s, where the expression for Im s Ti, ... 
was substituted. Note that one has to perform one subtraction in order to render the 
dispersion integral convergent. 

In addition, we applied a new method which completely avoids the use of explicit 
imaginary parts of hard-scattering amplitudes, allowing one to numerically calculate the 
NLO parts of LCSR, e.g., F\(q 2 , M 2 , Sq ) in ( [4.7D , analytically continuing integrals to the 
complex plane. We make use of the fact that the hard-scattering amplitudes Ti,Tf , Tf are 
analytical functions of the variable s = (p + q) 2 in the upper half of the complex plane, 
because of ie's in Feynman propagators. Consider, as an example the twist-2 part of F± 
given by the integral over s in the second line of ( [4.71 ). Since the integration is performed 
along the real axis, the operation of taking the imaginary part can be moved outside the 
integral. To proceed, one has to shift the lower limit of the s-integration to any point at 
q 2 < s < m 2 . This is legitimate because all Xi's are real at s < m 2 . Then one deforms 
the path of the s-integration, replacing it by a contour in the upper half of the complex 
plane, as shown schematically in Fig. [| so that all poles and cuts are away from the 
integration region. Only when s is approaching the upper limit Sq , one nears the pole 
at u = (m 2 — q 2 )/(sQ — q 2 ) while performing the integration over u. Because this pole 
does not touch the limits u = 0, 1, it is possible to avoid it by moving the contour of the 
it-integration into the upper half of the complex it-plane (see Fig. |3|). After that, both 
numerical integrations become completely stable. Note, that in both s- and u- integrations, 
we integrate over the semi-circle, but the contour of the integration can be deformed in 
an arbitrary way in the upper half of the complex plane. The numerical integrations of T\ 
over these contours yield an imaginary part which represents the desired answer for F\ . We 
have checked that the numerical results obtained by this alternative method coincide with 
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twist 


Parameter 


Value at /i = 1 GeV 


Source 


2 


a 2 
a >4 


(J. 25 ± U.15 

-af + (0.1 ±0.1) 

n 
u 


average from pi 
7T77* form factor 


\ 

30] 


3 


/3tt 


1.74ig;£ GeV 
0.0045 ± 0.0015 GeV 2 
-1.5 ±0.7 


GMOR relation; m u ^ f 
2-point QCD SR 
2-point QCD SR[ 


ro 
19 
19 


m H 

1 


4 


« 


0.18 ±0.06 GeV 2 
f (0.2 ±0.1) 


2-point QCD SR 
2-point QCD SR 


19 
19 


1 
1 



Table 1: Input parameters for the pion DA 's. 



the ones obtained by the direct integration over the imaginary parts, thereby providing an 
independent check. 

5. Numerical results 



Let us specify the input parameters entering the LCSR fl4.4D , ( |4.&j ) and ( 4.1l| ) for B — > tt 
form factors and the two-point sum rule (|C.1| ) for 

The value of the 6-quark mass is taken from one of the most recent determinations 



m b (m b ) = 4.164 ± 0.025 GeV , (5.1) 



based on the bottomonium sum rules in the four-loop approximation. Note that (|5.1| ) has 
a smaller uncertainty than the average over the non-lattice determinations given in p8[ : 
mb(mb) = 4.20 ± 0.07 GeV. However, as we shall see below, the uncertainty of m^mi,) 
does not significantly influence the "error budget" of the final prediction. Furthermore, 
in our calculation, the scale-dependence m;,(// m ) is taken into account in the one-loop 
approximation which is sufficient for the 0(a s )-accuracy of the correlation function. Note 
that using the MS mass inevitably introduces some scale-dependence of the lower threshold 
m 2 in the dispersion integrals in both LCSR and /b sum rule. However, this does not 
create a problem, because the imaginary part of the OPE correlation function obtained 
from a fixed-order perturbative QCD calculation is not an observable, but only serves as 
an approximation for the hadronic spectral density. 



The QCD coupling a s (fi r ) is obtained from a s (mz) = 0.1176 ± 0.002 [28|, with the 
NLO evolution to the renormalization scale fi r . In addition to [i m and /i r , one encounters 
the factorization scale [if in the correlation function, at which the pion DA's are taken. In 
what follows, we adopt a single scale fi = ii m = fJ> r = in both LCSR and two-point SR 
for fs- The numerical value of [i will be specified below. 



The input parameters of the twist-2 pion DA include = 130.7 MeV [28] and the 
two first Gegenbauer moments a 2 and aj normalized at a low scale 1 GeV. For the latter 
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we adopt the intervals presented in Table [p. The range for 0^(1 GeV) is an average [19] 
over various recent determinations, including, e.g., a^lGeV) = 0. 26+009 calculated from 
the two-point sum rule in pjJ. For a\ we use, following [pX] ], the constraint a^lGeV) + 
aJ(lGeV) = 0.1 ± 0.1, obtained from the analysis of ^77* form factor. Having in 
mind, that at large scales the renormalization suppresses all higher Gegenbauer moments, 
we set a> 4 = in our ansatz for ^^(u) specified in App. A. The uncertainties of aJ^lGeV) 
remain large, hence we neglect very small effects of their NLO evolution taken into account 
in |. 

The normalization parameter fj, n (lGeV) of the twist-3 two-particle DA's presented in 
Table [l] is obtained adopting the (non-lattice) intervals [^] for the light quark masses: 
m u (2 GeV) = 3.0 ± 1.0 MeV, m d (2 GeV) = 6.0 ± 1.5 MeV. Correspondingly, the quark- 
condensate density given by GMOR relation is: 

<g<?)(lGeV) = -i/^(lGeV) = -(246±?| MeV) 3 , (5.2) 

where very small 0(m^ d ) corrections are neglected. We prefer to use the above range, 
rather than a narrower "standard" interval (qq) (lGeV) = -(240 ± 10 MeV) 3 employed 



in the previous analyses. In fact, (5.2) is consistent with (qq) (lGeV) = (254 ± 8 MeV) 3 
quoted in the review pi|| , as well as with the recent determination of the light-quark 
masses from QCD sum rules with 0(a 4 ) accuracy |H§: m u (2 GeV) = 2.7 ± 0.4 MeV, 
m d (2 GeV) = 4.8 ± 0.5 MeV. 

The remaining parameters of the twist-3 DA's (f^, uj^-k) and twist-4 DA's (5%, %) 
presented in Table || are taken from [19], where they are calculated from auxiliary two- 



point sum rules. The latter are obtained from the vacuum correlation functions containing 
the local quark-gluon operators that enter the matrix elements (A.6), ( |A.7| ) and ( |A.12j ), 



( A.13j ). The one-loop running for all parameters of DA's is taken into account using the 



scale-dependence relations presented in App. A. Note that the small value of effec- 
tively suppresses all nonasymptotic and three-particle contributions of the twist-3 DA's. 
Furthermore, the overall size of the twist-4 contributions to LCSR is very small. Hence, 
although the parameters of the twist-3, 4 DA's have large uncertainties, only the accuracy 



of plays a role in LCSR 1 . Finally, in the sum rule ( C.l ) for fg the gluon condensate den- 
sity (a s /irGG) = 0.012^q'q^2 GeV 4 and the ratio of the quark-gluon and quark-condensate 
densities m 2 , = 0.8 ± 0.2 GeV 2 |H|] are used, the accuracy of these parameters playing a 
minor role. 

The universal parameters listed above determine the "external" input for sum rules. 
The next step is to specify appropriate intervals for the "internal" parameters: the scale 
H, the Borel parameters M and M and the effective thresholds s$ and Sq . In doing that, 
we take all external input parameters at their central values, allowing only and a\ to 
vary within the intervals given in Table |[ 

From previous studies |, |, |, |HJ it is known that an optimal renormalization scale 
is ~ \ m 2 B — m 2 ~ \J2rn^K (where A does not scale with the heavy quark mass), 



1 We also expect that the use of the recently developed renormalon model ]3^] for the twist-4 DA's, 
instead of the "conventional" twist-4 DA's 134] used here, will not noticeably change the numerical results. 
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and simultaneously, \x has the order of magnitude of the Borel scales defining the average 
virtuality in the correlation functions. In practice, M and M are varied within the "working 
windows" of the respective sum rules, hence one expects that also \x has to be taken in a 
certain interval. 

Calculating the total Borel-transformed correlation function (that is, the Sq — > oo 
limit of LCSR) we demand that the contribution of subleading twist-4 terms remains very 
small, < 3% of the LO twist-2 term, thereby diminishing the contributions of the higher 
twists, that are not taken into account in the OPE. This condition puts a lower bound 
M 2 > M^ in = 15 GeV 2 . In addition, in order to keep the a s -expansion in the Borel- 
transformed correlation function under control, both NLO twist-2 and twist-3 terms are 
kept < 30% of their LO counterparts, yielding a lower limit fi > 2.5 GeV. Hereafter a 
"default" value /i = 3 GeV is used. 

Furthermore, we determine the effective threshold parameter s B in LCSR for each 
M 2 > M^ in . We refrain from using equal threshold parameters in LCSR and two-point 
sum rule for fs, as it was done earlier, e.g. in @, ||. Instead, we control the duality 
approximation by calculating certain observables directly from LCSR and fitting them to 
their measured values. Importantly, we include in the fitting procedure not only Sq , but 
also the two least restricted external parameters and a\ , under the condition that both 
Gegenbauer moments remain within the intervals of their direct determination given in 
Table § 

The first observable used in this analysis is the S-meson mass. In a similar way, as 
e.g., in [10, [Pj| l, m 2 B is calculated taking the derivative of LCSR over — 1/M 2 and dividing it 



by the original sum rule. The .B-meson mass extracted from LCSR has to deviate from its 
experimental value mg = 5.279 GeV by less than 1 %. Secondly, we make use of the recent 
rather accurate measurement of the ^-distribution in B — > ttIv by BABAR collaboration 
1 35]. We remind that LCSR for B — > it form factors are valid up to momentum transfers 



q 2 ~ m 2 - 2mfeA , typically at < q 2 < 14 - 15 GeV 2 . To be on the safe side, we 
take the maximal allowed q 2 slightly lower than in the previous analyses and calculate 
the slope /^ 7r (g 2 )// i ^ 7r (0) from LCSR at < q 2 < 12 GeV 2 . The obtained ratio is then 
fitted to the slope of the form factor inferred from the data. We employ the result of 



1 36], where various parameterizations of the form factor fg~(c[) are fitted to the measured 



(^-distribution. Since all fits turn out to be almost equally good, we adopt the simplest 



BK-parameterization [37]: 



f +(BK) (a 2 ) 1 

{q ' - (5.3) 



f+ (BK \0) (1 " 9 2 /m|0(l " a BK q 2 /m 2 B ) 



with the slope parameter om = 0.53 ± 0.06 from |36| (close to olbk fitted in [p5|[). 

After fixing Sq for each accessible M 2 , we demand that heavier hadronic states con- 
tribute less than 30% of the ground-state B meson contribution to LCSR. This condition 
yields an upper limit M 2 < M^ ax = 21 GeV 2 . The resulting spread of the threshold 
parameter and Gegenbauer moments when M 2 varies between M^ in and M^ ax is very 
small: sfi = 36 - 35.5 GeV 2 , a^(lGeV) = 0.15 - 0.17, aJ(lGeV) = 0.05 - 0.03. The 
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f Bjl + (g 2 ) /f BjT + (0) 1 - 
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Figure 4: The LCSR prediction for the form factor shape /ij 7r (9 2 )//B 7r (0) fitted to the BK param- 
eterization of the measured q 2 -distribution. The two (almost indistinguishable) curves are the fit 
and the parameterization (5.5) at olbk = 0.53. 



quality of the fit is illustrated in Fig. || where the two curves: the calculated g 2 -shape of 



the form factor and the BK-parameterization (5J3) are almost indistinguishable. Thus, in 
our numerical analysis we "trade" the c/ 2 -dependence predicted from LCSR for a smaller 
uncertainty of the Gegenbauer moments and for a better control over the quark-hadron 
duality approximation. 

In the final stage of the numerical analysis we turn to the two-point QCD sum rule for 
Jb presented in App. C and find that at the adopted value of the renormalization scale [i = 3 
GeV the interval M = 5.0± 1.0 GeV 2 satisfies the same criteria as the ones imposed in the 
numerical analysis of LCSR: the smallness of higher power terms in OPE, and suppression 

^-0.9 
'+2.1 



of the heavier hadronic contributions. The threshold parameter Sq — ,,; " ' f-" vJ 



fixed by calculating m B from this sum rule. This time the deviation from the experimental 



value is even less than 0.5 %. For completeness, we quote the resulting interval fg = 214 , 5 



MeV. Note that the 0(a 2 ) correction taken into account in [18] is not included here. As 
usual, employing the sum rule for fg in order to extract the form factor from the LCSR for 
the product /s/g turns out to be extremely useful. One observes a partial cancellation 
of the a s -corrections in both LCSR and two-point sum rule and a better stability with 
respect to the variation of scales. 

To demonstrate some important numerical features of the LCSR prediction, in Fig. || 
(left) we plot the M 2 -dependence of the form factor fg~(0) with all other inputs fixed at 
their central values. The observed stability, far beyond the adopted "working" interval in 
M 2 , serves as a usual criterion of reliability in QCD sum rule approach. The ^-dependence 
plotted in Fig. || (right) is very mild from \i = 2.5 GeV up to \i = 6 GeV. The numerical 
size of the gluon radiative corrections in LCSR is illustrated in Fig. ||. 

The numerical analysis yields the following prediction for the vector B — > it form factor 
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Figure 5: Dependence of fg (Q) on the Borel parameter (left) and renormalization scale (right). 
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Figure 6: Gluon radiative corrections to the twist-2 (dotted line) and twist-3 (solid line) parts of 
LCSR for fg^iq 2 ), as a function of q 2 . The part proportional to (/) a (<p p ) is shown separately by 
dashed (dash-dotted) line. 



at zero momentum transfer: 



/+J0) = 0.263 



-0.004 
-0.005 



M,M 



+0.009 
-0.004 



±0.02 



shape 



+0.03 
-0.02 



±0.001 



(5.4) 



nib 



where the central value is calculated at fi = 3.0 GeV, M 2 = 18.0 GeV 2 , Sq = 35.75 
GeV 2 , a^(lGeV) = 0.16, aJ(lGeV) = 0.04, M 2 = 5.0 GeV 2 and sfi = 35.6 GeV 2 . The 
percentages of different contributions to the central value in ( |5.4[ ) are presented in Table |2[ 
In ( [5.4D the first (second) uncertainties are due to the variation of the Borel parameters 
M and M (scale /i) within the intervals specified above. The third uncertainty reflects the 
error of the experimental slope parameter. In addition, we quote the uncertainties due 
to limited knowledge of the "external" input parameters. We have estimated them by 
simply varying these parameters one by one within their intervals and fixing the central 
values for all "internal" input parameters. Interestingly, the largest uncertainty of order 
of 10% is due to the error in the determination of light-quark masses transformed into 
the uncertainty of the coefficient of the large twist-3 LO contribution. The spread 
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caused by the current uncertainty of rrib(mb) is much smaller, hence, does not influence the 
resulting total uncertainty, even if one increases the error of m^-determination by a factor 
of two. Remaining theoretical errors caused by the current uncertainties of a s , twist-3,4 
DA's parameters and higher-dimensional condensates are very small, and for brevity they 
are not shown in (5.4). 

Finally, we add all uncertainties in 
quadrature and obtain the interval: 



6-quark mass 


MS 


pole 


input 


central 


set II from pfl 


J&(0) 


0.263 


0.258 


tw2 LO 


50.5% 


39.7% 


tw2 NLO 


7.4% 


17.2 % 


tw3 LO 


46.7% 


41.5 % 


tw3 NLO 


-4.4% 


2.4 % 


tw4 LO 


-0.2% 


-0.9% 



26 +0 - 04 



(5.5) 



Table 2: The form factor f^ at zero momen- 
tum transfer calculated from LCSR in two different 
quark-mass schemes and separate contributions to 
the sum rule in %. 



which is our main numerical result. It can 
be used to normalize the experimentally 
measured shape, e.g., the one in (^), 
yielding the form factor f^ w (q 2 ) in the 
whole g 2 -range of B — > irlvi. 

With this prediction at hand, we are 
in a position to extract \V u b\. For that we 
use the interval 



9.1 ±0.61 



shape 



±0.3 



BR 



x 10" 



inferred [36| from the measured g 2 -shape |3 
|3||. We obtain: 



and average branching fraction of B 



\V ub \ 



3.5 ±0.4 ,. ±0.2 , ±0.1 L 

lift I shape \ BR 



x 10" 



(5.6) 
(5.7) 



where the first error is due to the estimated uncertainty of fg n (0) in ( |5.5[ ), and the two re- 
maining errors originate from the experimental errors in ( |5.6| ) . A possible small correlation 
between the shape uncertainty of our prediction for the form factor and the experimental 
shape uncertainty is not taken into account. 

The remaining two B — > ir form factors can now be predicted without any additional 
input. In particular, we adopt the same Borel parameter M 2 and effective threshold Sq , 
assuming that they only depend on the quantum numbers of the interpolating current for 
B meson. The scalar form factor f^ w (q 2 ), obtained by combining the LCSR for fi and 
(/btt+Zbtt)) ano - ^ ne penguin form factor f^ n {q 2 ) are presented in Fig. 0, in comparison with 
f^{q 2 )- The predicted interval for the penguin form factor at zero momentum transfer is 



/^(O) = 0.255 ±0.035, 
adopting /i = 3 GeV as the renormalization scale of the penguin current. 



(5.8) 



6. Discussion 



In this paper, we returned to the LCSR for the B — > it form factors. We recalculated the 
0{a s ) gluon radiative corrections to the twist-2 and twist-3 hard-scattering amplitudes and 
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Figure 7: The LCSR prediction for form factors f^(q 2 ) (solid line), f%^{q 2 ) (dashed line) and 
fs^iq 2 ) (dash-dotted line) at < q 2 < 12 GeV 2 and for the central values of all input parameters. 



presented the first complete set of expressions for these amplitudes and their imaginary 
parts. For the radiative corrections to the twist- 2 part of the LCSR for fg we reproduced 
the results of [||, ||. For the radiative corrections to the twist-3 part we confirmed the 
cancellation of infrared divergences observed in || in the case of asymptotic DA's. 
Including the nonasymptotic effects in these radiative corrections demands taking into 
account the mixing between two- and three-particle DA's. In fact, the parameter 
determining the size of nonasymptotic twist-3 corrections is numerically small, hence these 
corrections are not expected to influence the numerical results. 



Throughout our calculation and in the final sum rule relations we used the MS'-mass 
of the b quark, which is the most suitable mass definition for short-distance hard-scattering 
amplitudes. Indeed, as follows from our numerical analysis, the 0(a s ) corrections to the 
sum rules turn out to be comparably small. To demonstrate that, we returned to the 
pole-mass scheme in LCSR and used exactly the same input as in [[[(]] (the preferred "set 
2" with m^ ole = 4.8 GeV). We calculated the total form factor and separate contributions 
to the sum rule in both quark-mass schemes and compared them in Table 2. Note that 
the twist- 2 NLO correction is distinctively smaller in the MS'-mass scheme. In the twist-3 
part of the sum rule the a s -correction is small in both schemes. In the MS scheme, as seen 
from Fig. ^ this correction is dominated by the contribution of the DA ^ . In the pole 
scheme there is a partial cancellation between the contributions of the two twist-3 DA's. 
Our numerical results for the form factors f^ ,T (q 2 ) in the pole scheme are very close to 
the ones obtained in [10]. It is however difficult to compare separate contributions, because 



they are not presented in [10|. We also cannot confirm the numerical values of the twist-3 
NLO corrections to the form factor f^ n (q 2 ) plotted in the figure presented in the earlier 
publication Q. 

Further improvements of LCSR are possible but demand substantial calculational ef- 
forts. For example, obtaining radiative corrections to the three-particle twist-3, 4 contribu- 
tions is technically very difficult. Again, we expect no visible change of the predicted form 
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factors because three-particle terms are already very small in LO. A more feasible task is 
to go beyond twist-4 in OPE and estimate the twist-5,6 effects, related to the four-particle 
pion DA's, at least in the factorization approximation, where one light quark-antiquark 
pair is replaced by the quark condensate (In LCSR for the pion electromagnetic form 
factor these estimates have been done in ]39[). 

The numerical analysis of LCSR was improved due to the use of the g 2 -shape measure- 
ment in B — > ixlv. A smaller theoretical uncertainty of LCSR predictions can be anticipated 
with additional data on this shape, as well as with more accurate determinations of b- and, 
especially, u, ci-quark masses. 

In this paper, all calculations have been done in full QCD with a finite 6-quark mass. At 
the same time, the whole approach naturally relies on the fact that is a very large scale 
as compared with Aqcd and related nonperturbative parameters. Our results demonstrate 
that the twist-hierarchy as well as the perturbative expansion of the correlation function 
work reasonably well. An interesting problem is the investigation of the — > oo limit of 
LCSR and various aspects of this limiting transition, e.g., the hierarchy of radiative and 
nonasymptotic corrections. This problem remaining out of our scope was already discussed 
in several papers: earlier, in |7| at the LO level, in ||,[40] at NLO level and more recently, 



in ill 



[ref.] 


fg n {q 2 ) calculation 


IbM 2 ) in P ut 


\V ub \x 10 3 




lattice (rif = 3) 




3.78±0.25±0.52 


El 


lattice (rif = 3) 




3.55±0.25±0.50 


[43] 




lattice SCET B -> vrvr 


3.54 ±0.17 ±0.44 


[44] 




lattice 


3.7 ±0.2 ±0.1 


\E 




lattice LCSR 


3.47 ±0.29 ±0.03 


[10, 36 


LCSR 




3.5 ±0.4 ±0.1 


this work 


LCSR 




3.5 ±0.4 ±0.2 ±0.1 



Table 3: Recent \ V u b\ determinations from B — > irlvi 



Finally, in Table ||| we compare our result for \V u b\ with the one of the previous LCSR 
analysis and with the recent lattice QCD determinations obtained at large q 2 and ex- 
trapolated to small q 2 with the help of various parameterizations. The observed mutual 
agreement ensures confidence in the continuously improving V u b determination from exclu- 
sive B decays. 
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A. Pion distribution amplitudes 



For convenience, we specify the set of the pion DA's and their parameters used in this 
paper. The notations and parameters for twist-3 and 4 DA's are taken from [pj, where 
the earlier studies [34, are updated. 

The two-particle DA's of the pion enter the following decomposition of the bilocal 
vacuum-pion matrix element (for definiteness, ir + in the final state): 

(n+(p)\ui( Xl )4(x 2 )\0) x 2^ = jf * du e <«P-*i-NflP-*a ^ 75 ]^» 

1 " \ 

+-^\jhs]iiu{x 1 - x 2 ) 2 4>^(u) - - ^2)75]^ y ip4K(v)dv j , (A.l) 

' 

In the above, the product of the quark fields is expanded near the light-cone, that is, 
x% = £iX, where £j are arbitrary numbers, and x 2 = 0; u = 1 — u. The path-ordered 
gauge-factor (Wilson line) is omitted assuming the fixed-point gauge for the gluons. The 
light-cone expansion includes the twist- 2 DA <p n , two twist-3 DA's 0g , 4>^ n and two twist- 
4 DA's 4>4 n and ip^. The usual definitions of DA's are easily obtained, multiplying both 
parts of (A.l) by the corresponding combinations of 7 matrices and taking Dirac and color 
traces. 

The decomposition of the three-particle quark-antiquark-gluon matrix element is: 
(K + ( P )K(xi)g s G« u (x 3 )di(x 2 )\0) x ^ = ^ J p a ^K-i+^+«3*3) 

iP^X -Pn9v\)^4m{(y.i) H ^(p ■ x ) " ~ ( $47r ( Qi ) + ^(ai)) ) 



- l -y^v6 P (lx)^{(p P g SX -pV A )^(«i) + PX(PS *[ (<M«*) + * 47 r(0i)) } 



(A.2) 



including one twist-3 DA and four twist-4 DA's : $4^, ^4^, ^a-k and ^4^. Here the 
convention e 0123 = —1 is used, which corresponds to Tr{^^^^ u ^ a ^} = 4ie^ m ^ . 

The following expressions for the DA's entering the decompositions ( |A.1| ) and ( |A.2[) 
are used: 

• twist- 2 DA: 



(p n (u) = 6uu[l + a 2 C 2 (u — u) + a4C 4 (u — u)j, (A.3) 

where, according to our choice, the first two Gegenbauer polynomials are included in 
the nonasymptotic part, with the coefficients having the following LO scale depen- 
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dence: 

25C F 91C F 

02(^2) = [L(P2, Mi)] 6/3 ° fl^Mi), 04(^2) = [£(M2, A*i)] 15/30 04 (/ii) (A.4) 
with L(/x 2 ,Mi) = a s (fj, 2 )/a s (fj,i), f3 = 11 - 2n//3. 
• twist-3 DA's : 



^37r(ai) = 360aia 2 ai 



l + ^(7a 3 -3) (A.5) 



2 

with the nonperturbative parameters ^3^ and defined via matrix elements of the 
following local operators: 



{n (p)\u<7nvl5G a pd\0) = if37r[(PaP»g/3u ~ Pj3Pix9au) ~ (PaPu9p^ ~ P/3Pv9a»)\ ,(A.6) 

3 3 
(7r + (p)|ucr^ A 7 5 [Z) /3 , G aX ]d - -dpua^xlbGaxdlO) = -—h^^PaPpP^. • (A.7) 

The scale dependence of the twist-3 parameters is given by: 



1 1 'v F 



M^) = [L(fi 2 ,A*l)] ^M^l); /a* 0*2) = [L(n2,fJ-l)] fio ^ 3 ' /37r(Ml) ( A - 8 ) 
(/ 3 ^)(^2) = [L(^2,/"l)]^^ +10 ) (/3*W!br)(Ml) . (A.9) 

The corresponding expressions for the twist-3 quark-antiquark DA are: 
= 1 + 30^C 2 1/2 ( U -a)- 3^^C 4 1/2 (n - «), 
O) = 6u(l - «) (l + 5^- (l - ^p) C 2 3/2 (n - «)) . (A.10) 
• twist-4 DA's: 

$47r(ai) = I205le w (ai - a 2 )aia 2 a 3 , 

^4n(ai) = 305l(fi){ ai - a 2 )a 2 3 [^ + 2e n (l - 2q 3 )] , 

$47r(ai) = -1205laia 2 a 3 [^ + e n (l - 3a 3 )] , 

* 47r (ai) = 30<5 2 a|(l - a 3 )[^ + 2^(1 - 2a 3 )] , (A.11) 

are the four three-particle DA's, where the nonperturbative parameters 8 2 and e n are 
defined as 

(7r + (p)\uG a ^ a d\0) = iftUpp , (A.12) 
and (up to twist 5 corrections): 

(n + (p)\u[D^ G u ^d - ^uG ua Zd\0) = ~U^p^p v , (A.13) 
with the scale-dependence: 

8C F 10 

S 2 M = [L(//2, A*i)] 35 ° ^(/^i) , (fen)(j*2) = [L(a*2,Mi)]*> (<£<*)(/*i) • (A.14) 
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Note that the twist-4 parameter lo^ introduced in [19| is replaced by = (21/8)tJ4 7r . 
Correspondingly, the two-particle DA's of twist 4 are: 

47r (it) = — 5lu 2 u 2 + S5leJuu{2 + 13mt) + 2u 3 (10 - Ibu + 6u 2 ) kiu 

+2n 3 (10- 15u + 6n 2 )lnu} , (A.15) 

^M«) = y^C 2 1/2 (2n-l). (A.16) 



These DA's are related to the original definitions 34] as 

u 

dg2{u) 



<M«) = 16(<7i(u) - / g2(v)dv), ^(u) = -2^^ • (A.17) 



B. Formulae for gluon radiative corrections 

Here we collect the expressions for the hard-scattering amplitudes entering the factoriza- 
tion formulae ([OJ) and the resulting imaginary parts of these amplitudes determining the 
radiative correction ( [Q| ) to LCSR ( fO| ) for fi, as well as the analogous expressions for 
LCSR and ( 4.1 1| ) for the other two form factors. 

To compactify the formulae, we use the dimensionless variables 

ri = 4, r2 = &±|£, (B.l) 
(in the imaginary parts r 2 = s/m 2 ) and the integration variable : 

1 r 2 



p = n + u(r2 — n) [ du = f — , (B.2) 

J J r 2 -n 



''1 



and introduce the combinations of logarithmic functions 

G(x) = Li 2 (x) + ln 2 (l -x)+ ln(l - x) fin - 1 V (B.3) 
where Li2(x) = — T ~~ * s ^ ne Spence function, and 

L.W-tafCjZ^)-!, i , W = lll ffez>£4)_I. (B.4) 
\ x p z J \ x n £ J x 

The imaginary parts of the hard-scattering amplitudes are taken at fixed g 2 < m 2 
(n < 1), analytically continuing these amplitudes in the variable s = (p + q) 2 (or r 2 ). The 
result contains combinations of 8(1 — p), 6(p — l) and S(p— 1) and its derivatives. To isolate 
the spurious infrared divergences which one encounters by taking the imaginary part, we 
follow M and introduce the usual plus-prescription 
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for generic functions (p(p), <?(p). Furthermore, to make the formulae for imaginary parts 
more explicit, we partially integrate the derivatives of 8{p — 1) using, e.g.: 



jT 2 dp5'(p - 1)<Xp) = J d P 5(p - 1) (--jL + 6(r 2 - 1)) (f>(p) , (B.6) 



omitting the terms with 5(r 2 — 1) in all cases where (f)(1) = 0. 
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As already mentioned, the above formulae are obtained in the MS scheme. To switch 
to the one-loop pole mass of b quark the following expressions 
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have to be added to the hard-scattering amplitudes Ti, Tf and Tf . respectively. The 
corresponding additions to the imaginary parts are 
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B.2 Amplitudes for + /~J LCSR 
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B.3 Amplitudes for f% n LCSR 
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C. Two-point sum rule for fs 



(B.39) 
. (B.40) 



We use the sum rule with 0{a s ) accuracy with the perturbative part calculated in the MS 
scheme for 6-quark [O] : 
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(C.l) 



where M and are, respectively, the Borel parameter and effective threshold. In the 
above, the functions determining the spectral density of the 0(a s ) corrections to the per- 
turbative and quark condensate terms are 



Pl (s,m 2 b ) = ^{l-x)l (1 -aj)[4Li 2 (a;) + 21nxln(l - x) - (5 - 2x)ln(l - x)] 



+(1 - 2x)(3 - x) Inx + 3(1 - 3x) In ( ) + ^(17 - 33x) } , (C.2) 



mr 



where x = — , and 



<5i(M 2 ,m 2 ) 



AT/ V M 



hi 



m 2 ^ 3 



(C.3) 



respectively, and T(n, z) is the incomplete T function. 
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